The optical conductivity of half-filled Hubbard ladders 
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We investigate the optical conductivity of half-filled N-leg Hubbard ladders far into the "decon- 
finement" limit (i.e., weak Hubbard interaction and relatively strong interchain hopping). The N-leg 
Hubbard ladder is equivalent to an N-band model with velocities obeying vi = vn < V2 — Wjv-i < .... 
When N is not too large (N = 3, 4...), the band pairs (i,N + 1 — i) successively flow to the D-Mott 
state leading to a cascade of charge and spin gaps [U. Ledermann, K. Le Hur, and T. M. Rice, 
Phys. Rev. B 62, 16383 (2000)], and to the progressive closing of the two-dimensional (2D) Fermi 
surface (FS). The optical conductivity at finite temperatures can then exhibit coexistence between 
a prominent Drude peak and a high-frequency preformed pair continuum, split by sharp excitonic 
peaks arising due to an approximate SO(8) symmetry. For very large (but finite) N, all neighboring 
bands interact on the 2D FS, leading to a low-temperature 2D Mott crossover accompanied by a 
Spin Density Wave (SDW) instability (similar to the 2D case) . In this limit the optical conductivity 
exhibits a unique charge gap -excitations above the gap being bound hole-pairs- and that the exciton 
features vanish. These results could help to explain the optical conductivity of 2D systems at and 
close to half-filling, an example of which is the pseudogap phase of high-T c cuprates. 

PACS numbers: 71.10.Pm;71.30.+h;72.10.-d 



I. INTRODUCTION 

Understanding the nature of the interplay between 
strong interactions (generally emerging at the stoichio- 
metric electron density due to the possibility of umk- 
lapp scattering) and dimensionality is an important open 
matter relevant to a large class of materials including 
the cuprate superconductors and quasi one-dimensional 
(ID) organic conductors. The latter can be viewed as 
possible realizations of coupled one-dimensional Hubbard 
chains at half-filling where electrons can hop from chain 
to chain (although the quarter-filled chains are weakly 
dimerized allowing half-filled umklapps, there is increas- 
ing evidence that quarter-filled umklapp scattering dom- 
inates the physics) 1 . In particular, the Fabre salt family 
(TMTTF) 2 X (X = C10 4 , Br and PF 6 ) displays insulat- 
ing behavior at ambient pressure up to quite high tem- 
peratures. This reflects the presence of strong umklapp 
scattering (interaction, U) along the chains due to the 
commensurate filling resulting in a ID Mott transition; 
essentially one has virtually uncoupled insulating chains 
(confinement 

By contrast, substitution of Se for S to create the Bech- 
gaard salt (TMTSF) 2 PF 6 , increases the hopping between 
chains (tj_) or alternately decreases the dimcrization suf- 
ficiently to delocalize particles in the transverse direc- 
tion (inducing deconflnement). The system thus ex- 
hibits a crossover to a regime of metallic planes (as tj_ 
scales to strong coupling faster than U). Experimentally 
this dimensional crossover takes place around 100K in 
(TMTSF) 2 PF 6 ; This is manifested by a change from T 
(single-chain Luttinger liquid) to T 2 (Fermi liquid) be- 
havior of the dc transport along the chain axis^. How- 
ever, such compounds are located quite close to the con- 
finement transition such that the optical conductivity of 
(TMTSF) 2 PFg exhibits many features in common with 



the single-chain Mott insulator^ (TMTTF) 2 PF 6 . Indeed, 
only 1% of the spectral weight 6 contributes to the dc 
(Fermi liquid) transport. Some theoretical attempts to 
reproduce this result (which bears a superficial resem- 
blance to our half-filled chains for small N (see Fig. 1)), 
have been performed recently using Dynamical Mean 
Field Theory (DMFT)i and RPA£. 

Here, we are additionally motivated to tackle the 
physics of 2D Cu-0 planes of cuprate materials at and 
close to half- filling where Mott physics ( "Mottness"— ) 
is also ubiquitous^; We systematically investigate the 
behavior of the optical conductivity for half-filled N-leg 
Hubbard ladders (=coupled chains) far into the "decon- 
flnement" regime (the bare transverse hopping ampli- 
tude is almost equal to the longitudinal hopping am- 
plitude). For weak on-site repulsion, we can make use 
of our previous band approach for half- filled N-leg Hub- 
bard laddersii. This allows us to properly take into ac- 
count the renormalization of the different coupling chan- 
nels (umklapp, Cooper,...) and also to provide a rigorous 
study of the ground state. The N-leg Hubbard ladder is 
equivalent to an N-band model where at half-filling the 
Fermi velocities obey v\ — vn < i> 2 = fjv-i < 

For small N (N = 3,4...), this model results in a cas- 
cade of energy scales^i, where band pairs (i, N + 1 — i) 
successively flow to the D-Mott state which possesses an 
enlarged SO (8) symmetry; This is reminiscent of the two- 
leg ladder behavior—. This leads to a complex behavior 
of the optical conductivity: a Drude peak and a high- 
frequency particle- hole continuum coexist, with sharp ex- 
citon peaks appearing below the continuum reflecting the 
underlying SO (8) symmetry. 

By increasing considerably the number of chains, 
"four-band" couplings (Fig. 2) become relevant on the 
2D Fermi surface. At sufficiently low temperatures, a 2D 
insulating transition with strongly enhanced SDW corre- 
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lations arises^. The low-energy physics already converges 
to that of the purely 2D Hubbard mode&&l*. We rigor- 
ously establish that charge excitations above the Mott 
gap consist of bound hole-pairs (preformed pairs) and 
that the excitonic peaks cannot survive in this limit. This 
could be relevant to an understanding of optical conduc- 
tivity measurements on the high-T c cuprates at and close 
to half-fillingi^; Indeed, umklapp scattering should not 
be ignored in the underdoped regime of the cuprates- 
features seen above ~ 2eV and the transfer of spectral 
weight as a function of doping or temperature (T) reflect 
this importance. 

II. THE MODEL 

Our starting point is the N-leg Hubbard ladder model, 
Hkitl — -t d] s (x + a)d is (x) + h.c. 

x,i,s 
x,i,s 

t and t± denote the hopping matrix elements along and 
perpendicular to the chains, a is a lattice step, and di S (x) 
annihilates an electron with spin s on chain i at the rung 
x. The interaction term reads 

Hi at = Uj2dU x ) d ir(^4i^)dii(x), (2) 

where U is the on-site Hubbard repulsion. Again, here 
we are investigating the perfect "deconfinement" regime 
which implies that the perpendicular hopping t±_ is suf- 
ficiently large to deconfine all the electrons in the trans- 
verse direction. More precisely, below, we consider the 
weak-interaction limit < U <C (t,t±) where tractable 
(and rigorous) calculations are indeed possible^. It is 
convenient to use the band picture where the kinetic part 
simply takes a diagonal form: 

H Km = [ dk e *( fc )*L( fc )*-( fc ); (3) 

»=1. ..«>•' 

and ^>i s are the creation and annihilation operators 
for the band i and 

= — 2t cos(fca) — 2t± cos(fcj_ia). (4) 

The transverse (Fermi) momenta here simply obey k±i — 
±wi/(a(N + 1)) and at half-filling the longitudinal mo- 
menta are exactly determined by Ci{kpi) = 0. The re- 
sulting Fermi velocities = ^ sm(kpia) then take the 
form 

Vi = Us = y \/t 2 -{t ± cos[7ri/(iV + l)]} 2 , (5) 

where t = N + 1 — i. The transformation from chain to 
band (with open boundary conditions) readsife 



For more details on the method at half- filling, see Ref.lllT 
Two different limits arise. When TV is not too large, 
careful investigation of the coupled Renormalization 
Group equations (RGEs) of the interacting N-band prob- 
lem yields a hierarchy of energy scales, Tj ~ f e -avih/(aU) 
(a is a constant parameter of the order of 1), Ti > T% > 
... > T r -for N even r = N/2 and for N odd r = (JV+l)/2, 
where band pairs (i,t) decouple from all other bands and 
freeze out. This decoupling, due to band pair umklapp 
terms satisfying (kpi + Hfi) — K /a, opens both a charge 
and spin gap. The effective single-particle gaps are ap- 
proximately given by Aj — Tj. More precisely, the cou- 
plings of the band pairs scale towards the two- leg 
ladder fixed point with SO(8) symmetry^. To compute 
the optical conductivity in this limit, it will be crucial 
to exploit this enlarged symmetry. When iV becomes 
large (but finite), our method is still rigorous provided 
the energy difference between two neighboring bands is 
much larger than the largest relevant energy gap in the 
problem, i.e., te~ l ^ u . This imposes the condition^ U < 
t/]nN,t±/lnN. Here, coupling between neighboring 
bands leads to a 2D-like SDW instability (crossover) - 
opening a Mott gap- at the temperature^: 

N 

EN/2 ~~H 

-,„,, i=l "i , (7) 

Since we are considering relatively large bare values of t± 
(i.e., perfect nesting on the 2D FS), the resulting charge 
gap A ~ T s dw is homogeneous on the 2D Fermi surface. 
When substantially decreasing t± (approaching the con- 
finement transition), we expect that the physics would 
have a larger dependence on the transverse momenta&ii 
k±i. Finally, ground-state charge properties still resem- 
ble those of the small-N (two-band) limit . 

III. SMALL N: MAPPING TO SO(8) 

For small N, the low-energy Hamiltonian is then the 
sum of N/2 [(N-l)/2 for N odd] two-leg ladder Hamilto- 
nians if^j corresponding to the band pairs (i,T) plus the 
Hamiltonian of a single chain for N odd H^L. 

Let us consider the case = 3 for simplicity. Here, 
the relevant couplings of bands 1 and 3 flow towards the 
universal ratios of the two-leg ladder, and therefore the 
low-energy Hamiltonian Hi 3 can be rewritten in terms 
of an SO (8) Gross-Neveu model (as shown previously by 
Lin, Balents, and Fisher for the two-band modeli^) 

Hi, 3 = [ dx S.-iv 1 J2vlr*d x y a -g(J2*lT v y a ) 2 \. 
J ^ 0=1 0=1 ' 

We have implicitly introduced operators $ R/La f° r right 
and left moving fermions and the Pauli matrices t 1 act 
on right and left sectors. Even though the four Dirac 
fermions above are related to the charge and spin degrees 
of the two bands, they cannot be simply identifed as the 
original (real) fermions on each band; for a complete map, 
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see Ref. ^3 The low-energy physics emerging from bands 
1 and 3 depends on a single coupling g (whose bare value 
is of the order of U). It is straightforward to show that 
the interacting part of Hi 3 both opens a charge and spin 
gap equal to 2A 1 = Ag < > ~ 2T X . The ground 

state corresponds to the "D-Mott state" : a Mott insula- 
tor having short-range pairing correlations with approxi- 
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FIG. 1: Optical conductivity of the three-band model for 
T 2 < T < Ti. The Fermi surface (FS) (drawn for t± = t) ex- 
hibits a blatant truncation, meaning that bands 1 and 3 form 
the D-Mott state opening a gap whereas the band 2 is still 
metallic (and described by a Luttinger theory). The optical 
conductivity then exhibits a Drude peak (due to band 2), a 
high-energy continuum (treated here at the mean field level) 
above the charge gap 2Ai and a sharp exciton peak due to 
the SO(8) symmetry of the D-Mott state. 

-mate d-wave symmetry 12 . Additionally fluctuations of 
the gap around its vacuum value can generate attrac- 
tive interactions between the fermions ty a leading to the 
formation of excitons, whose mass satisfies Mi = y/3Ai 
(i.e. is smaller than the charge gap of the (localized) 
preformed-pair continuum 12 ). Unlike the single-chainia, 
here excitons already appear for purely on-site repulsion. 

Band 2 is described by a single chain at half-filling 
(/1 — 0), which is embodied by gapless spinons (spin-1/2 
excitations) and the following charge Hamiltonian: 



m 



— {d x $> p f + K p (d x e p f 
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_9p_ 

(no) 



cos(v87r<I ) p) - nd x ®( 



(9) 



d x $ p describes charge density fluctuations and d x Q p cur- 
rent excitations. g p is of the order of U and u p is almost 
equal to v 2 . It is immediate to check that H 2 induces a 
charge gap 2A 2 ~ 2T 2 < 2Ti. Let us now compute the 
optical conductivity of the three-leg ladder. 

Using formula iJBJ, the electrical current can be written 
in the chain basis as: 



J(x) 



2i0F ■ 



2i^K 



^2jj-( d L( x + a)d ts {x) - d[ s (x)d is (x + a)) 



■Vdi 



(^l s (x + a)^ 2s (x) - ^l s (x)^2s(x + a)) 



+ VD1 ^ VD2 ^l(x + a)9i.{x) - + a) 

+^l{x + a)V 3s (x) - *l s (x)y 3s (x + a)) 

VDI - VD2 . Tff 



+ 



(^\ s (x + a)V 3s {x) + *l(x + a)* ls (x) 
-*\ s (x)>S> 3s (x + a) - 9l(x)* u {x + a))), (10) 



2n 



where we have summed up over the contribution coming 
from each spin s and each leg i = 1,2,3 of the ladder. 
Symmetry of the ladder dictates that the drift velocity 
(udi) in chains I and 3 be equal, although not necessarily 
the same as in chain 2 (vr>2), and the time interval for a 
hop of one lattice spacing, a, has been taken as 5ti = 
Note that were these velocities equal, the coefficients of 
the first two terms in the band basis would be identical 
while the third would disappear leaving a straight sum 
over bands. Indeed, the last term is found to contain 
an oscillatory contribution in x, which, when inserted 
into Eq. 12 vanishes. In linear response, the optical 
conductivity (at k = 0) takes the standard form 



5ie o{oj) = 3m 



n(w) 



to 



(11) 



the current-current correlator being 



1 



n(w) = -J dxdre^ T (T t J(x,t)J(0,0))^-«o+8- 
In the three-leg ladder, we find the following current 

\/2 sm(ltFia)d x Q p +) 



t ( a a , VD1 + Vm 
J = e[v D id x <d p H 

7T Z 



Jlband "f" J 2 bandi 



(13) 



where one can re-express the 2-band contribution in 
terms of Majorana fields (real chiral fermions ipR/L a — 

VR/L2a-1 + irjR/L2a) as 

J 2band = - {Vm t^ D2)e sHk F1 a)*iT^i 



2na 

(vdi +v D 2)ieV2 . , . 
= 1= sm(k F ia)(r] R iri R 2 -r]LiVL2)- (14) 

Using Eq. 12, the contribution of these terms to the 
conductivity is, 



, a , , -e 2 (v D i+v D2 ) 2 sm 2 (k F ia) 



mi 



2 hiu 



{ J dx 



dre WT <TtVpiVp2(x,t)7] P iVp2(0,0) > }. (15) 

P=R,L 



As pointed out in Ref. 12, since the two-band model is 
equivalent to the SO (8) Gross-Neveu model whose exci- 
tation spectrum is known, below the D-Mott gap mass 
Mi bound-states of charge ±2e fundamental fermions can 
form. The latter factor of Eq. 15 bears resemblance to 
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the (k = 0) Green's function of a localized particle of 
mass Mi, so that 



dx 



dTe iOx e iw, 



T <T TVl r) 2 (x,r)(mV2(0,0)) > 



dx 



dre Mx e lul " T Ae- 



« e 



A(t-£e-T) 



Mi 
n 



(16) 



where we have introduced a normalization constant, A, 
the length of the system, L, which we can take to be 
infinite relative to the correlation length, £ as we have 
power law correlations in this system. This means that 
for u> < 2Ai, the 2-band contribution to the optical con- 
ductivity at T=0 should be 



5tecr 2&and (w) = — 5(u3-Mx). 

a z nu> 

As in principle there is only one velocity scale in the (mas- 
sive) SO(8) Gross-Neveu model, v±, one might expect to 
find (vdi + VD2) <x v\. However, one might expect the 
opening of this gap to strongly renormalize v\ . 

An interesting situation emerges when the Fermi sur- 
face is partially truncated, i.e., when T2 <C T -C T\\ 
Band 2 is still metallic whereas bands 1 and 3 already 
form a D-Mott insulating state. In fact, the umklapp 



term 



9p 



j^jw cos(v87T$ p ) is still small, because its renor- 
malization has been (completely) cutoff by thermal ef- 
fects. Band 2 then behaves as a single-chain Luttinger 
liquid which produces a prominent Drude peak of height 
~ 2e 2 V£,i/h ~ 2e 2 u p K p /h at u> = 0; from this we 
see that Vdx — U2, the Fermi energy of the second 
band. Moreover, following Ref. 0, the current operator 
4'}r z, i'i clearly excites the mass V^Ai excitons, as well 
as higher energy continuum scattering states (preformed 
pairs) with energy above 2Ai. For clarity's sake, results 
have been summarized in Fig. 1. At very high-frequency, 
we approximately recover the behavior— of a single-chain 
at half-filling, i.e., <j{u>) cx u _1 (see Appendix A). 

IV. SMALL N: THE EFFECTS OF TEMPERATURE 

As mentioned previously, one certainly would expect a 
non-trivial role to be played by the temperature of the 
system. While at strictly T=0, one would expect to see 
a (^-function peak arising from the bound-state excitonic 
contribution to the optical conductivity, this peak should 
be thermally broadened. At T — 0, a mean field treat- 
ment would predict the spectral weight in the high-energy 
continuum to have a square-root singularity reflecting the 
van Hove singularity at the bottom of a band—:; At T 7^ 0, 
the scattering time between the fermions ^i becomes fi- 
nite and one would expect a(u> = 2Ai,T) » e Al ^ T (fi- 
nite). On the other hand, since the SO(8) Gross-Neveu 
model is integrable, the optical conductivity for this state 
can be found exactly in principle, as it consists of the 
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FIG. 2: Optical conductivity for (left) a 3-leg ladder, and 
(right) a 5-leg ladder. Here we have made the assumptions: 
Amid ~ lmeV w WK and t w 250meV, to yield f w 0.36, 
which defines: (left) Ai w 50.4^ and (right) Ai w 158^ 
and A2 w 20.9K. Further, temperature has been chosen to 
reflect the situation of minimal metallic conductivity: (left) 
T ss 207^ leading to approximately one third of the spectral 
weight in the Drude-like peak; and (right) T « 15K leading 
to approximately one fifth of the spectral weight in the Drude 
peak, and temperature broadening effects have been added in 
a minimal way (see text). 



bound-state, 2 particle scattering, and higher particle 
scattering contributions. Konik and Ludwig^S have cal- 
culated the (zero temperature) optical conductivity of the 
bound-state and the 2-particle scattering contribution of 
this model to be (k=0), 



exp 



dx G c (x) sinh 2 (x/3) 



+9(lu - Am 2 ) 



12m 2 



sinh(x) 



Am 2 



(to 2 — 3to 2 ) 2 u> 

x coshT 



■ exp 



dx 

x 



G c (x)(l — cosh(x) cos(- 



)) 



sinh(a;) 



]), (18) 



where A is an overall constant, m=Ai is the mass of 
an elementary SO(8) Gross-Neveu particle, and G c — 

-. In particular, it is useful to no- 



cosh(a;/2) 
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tice that the square-root singularity predicted by mean 
field theory does not survive in the exact calculation. In- 
tegration of their result, (assuming higher order scatter- 
ing processes do not dominate) shows that an astound- 
ing 69% of the spectral weight of the two-leg ladder (at 
zero temperature) is expected to lie within the bound- 
state. As noted, their results are exact up to the onset 
of the three-particle continuum, provided one stays at 
uj < lo c , where lu c is the temperature at which the RG 
couplings reach their approximate SO (8) coupling ratios. 
Above this frequency, presumably the chains will begin 
to behave more and more as individual chains, such that 
one might expect a crossover from the high frequency 
behavior found by Konig and Ludwig22 to an approxi- 
mate — as found by GiamarchiiS. The effects of such a 
crossover have been considered in Appendix A. 

In the hopes of making a more simple comparison to 
experimental curves, we here make the simplest assump- 
tions consistent with a low temperature broadening of the 
bound-state Gross-Neveu peaks following the treatment 
by Lin et ali&. There, it was argued that for temper- 
atures T<<m, one could essentially solve the classical 
Boltzmann equation to obtain a particle density due to 
thermal excitations of n ~ y/mTe^^ . It was further ar- 
gued that this approximation was self-consistent as the 
mean free path resulting I ~ i would be much greater 
than the de Broglie wavelength A ~ , 1 — , meaning 

° ° vmT 

that although the scattering processes themselves were 
quantum-mechanical, one could think of the particles es- 
sentially classically in between these time periods, lead- 
ing to a time between collisions t ~ fL . Furthermore, 
they argued that the shape of temperature-broadened 
bound-state peak should be proportional to this single- 
particle scattering time, additionally satisfying a scaling 
ansatz, <7 pea k ~ t£((c<j — y/3m)T) for (u> — v3m) << m. 

For our (N-leg) truncated ladders then, the simplest 
ansatz might be to assume that each <5-function peak 
is broadened by the lifetime of the particles associated 
with each successive gap on the Fermi surface; that 



is, to make the blanket assumption r. 



and 



represent the temperature broadened exciton peaks as 

In the case of interest, when 



"peak OC n 1 + ((w _y3 Ai)T .)2 

the temperature is above one of the gaps, one would ex- 
pect the spectral weight to be distributed about zero, 
but the breadth of the peak might still be governed by 
the ratio of the energy scales of the gap relative to the 
temperature. In such a way, we would expect to sim- 
ply shift the peak to zero, but effectively keep the same 
form as the excitonic peak broadenings. To create Fig. 
2, we have taken this ansatz in combination with the ex- 
act results of Konig and LudwigStt, and the assumption 
that as the ladder separates into effective 2- and 1-band 
ladder contributions, the spectral weight in each should 
be conserved. Additionally, we have input experimen- 
tal parameters in such a way as might be relevant for a 
discussion of (TMTSF) 2 PF 6 , although this leads to an 
unreasonably high value of H- w 0.36 for the validity of 



the SO (8) as noted in Appendix A. We note that these re- 
sults, even if one was well within the validity of the SO (8) 
symmetry, certainly underestimate the spectral weight in 
the continuum, as 3- and higher- particle processes have 
been neglected, in addition to the (unphysical) implicit 
assumption that there has not been a transfer of spec- 
tral weight between the bound-state and continuum as a 
function of temperature. 

None-the-less, the result is a semi-quantitative picture 
of what one might expect to find in the optical conduc- 
tivity of the ideal Hubbard ladder system at intermediate 
temperatures at half-filling (ie. if one could make mea- 
surements at frequencies sufficiently small frequencies 
and temperatures). The percentage of spectral weight 
in the Drude-peak would be expected to be diminished 
slightly as one would expect a small Mott contribution 
close to the relevant energy gap (here 1 meV), but the to- 
tal spectral weight of this feature when one has only one 
gap above the temperature scale would be expected to be 
of the order ■h , with a progressive narrowing of the width 
of this feature as the temperature decreased (as seen for 
example in (TMTSF) 2 PF 6 by Vescoli et a& For the 
small number of legs where this truncation should occur, 
one should not be able to achieve such a small percentage 
of the total weight (1%)& as seen experimentally^. 

V. LARGE N: ANTIFERRO MAGNETIC 4-BAND 
HOLE PAIRS 

Let us now increase considerably the number of chains 
(we consider the solvable case where the energy differ- 
ence between band pairs ~ 1/N is larger than the rel- 
evant gaps). As long as "4-band" interactions between 
bands (i, k,l,k) (fcj?, — > kFk,kF% — * k F j.) -which stand 
for 2D-like antiferromagnetic interactions (Fig. 3)- re- 
main small, our system is a bunch of (almost) decoupled 
bands and is still ID-like. However, at low-energy, the 
system clearly behaves as in 2D (See Ref. Indeed, 
two-band and single-band interactions will be dominantly 
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FIG. 3: Optical conductivity of the N-band model for large 
(but finite) N and T < T adw - 2D antiferromagnetic (AFM) 
interactions break the SO (8) symmetry (exciton peaks van- 
ish); The system flows to a 2D SDW state with a uniform 
charge gap. 
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renormalized by the "four-band" interactions: Certain 
couplings of the band pairs (i, i) still grow and tend to 
approach fixed ratios, however these differ from the small 
N ratios. Furthermore, for bands k and i which are close 
together on the FS, as k — > i, the four-band couplings 
become the same as the corresponding two-band cou- 
plings. The one-loop RGEs and theoretical details are 
given explicitly in Ref. Il3t To summarize, all the band 
pairs now strongly interact on the 2D FS, and this leads 
to a unique single-particle (Mott) gap A which can be 
calculated with the asymptotic ratios. Using Eq. J7J for 
t± — ► t, we find 

A~T sdw = te -*t/(U In N)_ (1Q) 

The logarithmic corrections come from the fact that v\ ps 
ta/N, a precursor of van Hove singularities in 2D. Let us 
repeat that here In AT < t/U for the (strict) validity of 
our calculations. In this limit, it is yet possible to derive 
the ground-state properties using bosonizationi^. 

The pinning of the antisymmetric spin mode between 
bands i and i, Q^a- ~ 0, produces spinon confinement 
leaving as physical particles spin-1 magnons. In the sym- 
metric spin sector, only the differences 8<jm+ — ® a kk+ 
are fixed, such that the total magnon mode(s), given by 

8s = #E^e ffS+ and $ s = ^/2/NY:?Ji$™+, 
remain gapless. The spin system behaves increasingly as 
a 2D SDW in agreement with the strong U limit 21 , and 
2D AFM interactions explicitly break the SO (8) symme- 
try of individual band pairs as it becomes more favor- 
able to pin the mode @ a „+ — CTfc fc+ rather than $ aii _^£. 
Finally, 2D AFM interactions do not much affect the 
ground-state charge properties of band pairs, i.e., as in 
the small-N limit $ P m+ « (symmetric charge mode) 
and Ops- ~ (antisymmetric superfluid mode). 

Let us now discuss the optical conductivity in this 
large-N limit. First, one can easily show that exciton 
peaks (that were a blatant signature of the underlying 
SO(8) symmetry for small-N) cannot persist. Integrat- 
ing out the spin sector -which now produces antiferro- 
magnetism and then totally decouples from the gapped 
charge sector at low energy- the symmetric charge mode 
Qpii+ is simply described by a Sine-Gordon modeU^ with 
j3 rs V4tt- The spectrum now contains only solitons and 
anti-solitons (bound hole- (electron) pairs) with disper- 
sion e(p) = ^Jp 1 + A 2 . At uj w and T < A, the 
density of excited carriers evolves like n ~ e~ A / T , which 
implies that the Drude peak has an exponentially vanish- 
ing weight, reflecting the simultaneous opening of a Mott 
gap on the whole 2D FS. Finally, at high-frequency, we 
still find o~(u>) oc lo^ 1 (we have almost decoupled bands). 

VI. CONCLUSIONS 

In closing, we are able to provide a microscopic (low 
energy) theory for half-filled N-leg Hubbard ladders - 
even when N is very large- taking into account both the 
ubiquity of Mott physics and the emergence of promi- 



nent antiferromagnetism. At small N, we have demon- 
strated that it is at least possible within the Hubbard 
model at half-filling for a system to exhibit metallic be- 
havior over a finite temperature range while maintaining 
a large Mott contribution, a situation not unlike that 
seen in (TMTSF) 2 PF 6 -although the T 2 resistivity en- 
countered in that case has not been recovered here, the 
relative spectral weight of the Drude feature would be 
O(jf) rather than 1%, and the relevant energy scales 
for this system are well outside the range of validity of 
our method. At high frequencies one should be able to 
employ a memory function approach to extract a de- 
tailed frequency dependence of the optical conductivity 
even in regions of moderate coupling. For large N (and 
U < tj In N, t±/ In N), the Mott gap opens simultane- 
ously on the 2D FS and, as a precursor of supercon- 
ductivity, charge excitations resemble hard-core bosons 
(preformed pairs), whereas spin excitations are gapless 
magnons (bosons). The optical conductivity is mostly 
sensitive to the Mott gap opening on the 2D FS (and not 
to the long-range magnetism) and excitons cannot per- 
sist! This seems to be in agreement with results on half- 
filled cupratesi^. When doping, we stress the occurrence 
of a spin gap and the progressive growth of d-wave phase 
coherence on the 2D FS. In order to properly discuss the 
pseudogap phase of cuprates, we must include the effect 
of a next nearest neighbor hopping 2 ^ t' . It is conceivable 
that such a term might provide a mechanism for the cre- 
ation of the "mid-gap" states of the doped cupratesi^S^ 
(which bear a qualitative resemblance to the two-leg lad- 
der contribution as outlined in Fig. 2), if the effective 
value of U were to vary substantially around the Fermi 
surface. 

We acknowledge discussions with C. Bourbonnais, R. 
Duprat, T. Giamarchi, S. Kancharla, A.-M. Tremblay, A. 
Tsvelik, and financial support by FQRNT and NSERC. 

APPENDIX A: GENERIC FEATURES FROM ONE 
TO THREE CHAINS 

A serious problem for this theory is the weakness of 
U required for the theory to flow to the strong coupling 
SO(8) fixed point. In particular, this imposes that our 
gap sizes and hence all of our temperature scales are out- 
rageously low, as the price to pay for the rigor of our 
treatment. In particular, to actually reach SO (8) within 
the validity of the 1 loop RG, the half-filled two- leg lad- 
der requires j < 2 x 1CU 6 , yielding a gap ~ £ e -500000^ 
which for all intents and purposes experimentally might 
as well be zero. This can be seen from Fig. 4, an explicit 
plot of the ratios of the couplings for the two-leg ladder 
as a function of the normalized number of steps. For a 
definition of the couplings of the two and three leg lad- 
der, the reader is referred to Refill and Ref.22. One sees 
that the curves described by the 1-loop RG equations are 
universal, although the region of validity decreases as U 
increases. In the inset, this is seen as a function of the 
maximal coupling, with ratios expected to reach 1 and 
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if the SO (8) fixed point is reached. 
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FIG. 4: Validity of the 1-loop RG: (main) Scaling of the RG 
equations as a function of the number of steps, normalized 
by the point at which approximate SO (8) symmetry is found 
(Imax). For a definition of the couplings of the two and three 
leg ladder, the reader is referred to Ref.ll and Ref.22. Cj 2 , 
fj 2 , and uj 122 have been plotted for U=10 -8 due to their 
close proximity, while the remaining curves have been plot- 
ted for t=o.09, u={f(r 8 ,i(r 6 ,icn 4 ,f(r 2 ,f} with \ max = 

(50993400, 508230, 5118.84, 51.7196, 0.51743) respectively. 
The vertical black dashed line delineates the point where the 
1-loop RG equations are expected to break down for U=10 -2 , 
the vertical gray dotted line for U=10 -4 . (inset) The same 

curves plotted as a function of the ratio of the coupling 
The vertical dashed line at 0.2, represents the point beyond 
which the 1-loop RG should not be valid. In principle, these 
curves could be used to extract the frequency dependence of 
the umklapp contributions to the two- (and three-) leg ladders 
as sketched in Fig. 7 and Fig. 8. 

Even in the absence of SO (8) symmetry, one clearly 
can say something about the evolution of the couplings 
(until the higher order contributions in U become rele- 
vant), which means that in principle one could calculate 
the w-dependence of the optical conductivity at high fre- 
quencies for the two and three leg ladders using a mem- 
ory function approach^ akin to Giamarchi's work on the 
single chain. 19 In this spirit, here we provide a heuristic 
sketch of such a derivation for the 3-leg ladder, and what 
one would expect to be the region of its applicability of 
this perturbative approach. At very high frequencies, one 
would expect the optical conductivity to break into two 
components roughly as, 
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,s = i "2e 2 u p K p | . *ti-u p+ j\ p+ . 2Q . 

H(U) + M lband (uj)) K(U) + M 2 band{u)) ' 

where the memory functions, Mu, an d(uj) and M 2 b an( z(u;) 
can be computed perturbatively in linear response theory 



as outlined in Ref. 19 and 24 respectively. In particular 
one finds that 19 in the limit u> 3> T, 
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For the three-leg ladder, the Luttinger coefficients can 



be written as, K„ — 
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lx , so that one sees that while K p 



In the high frequency limit, all 

7TI>1 |i + 

couplings, with the exception of uj 331 which vanishes so 
that one recovers the form of Ref. 24, are of order Ua/?l 
(if one is interested in the exact values, the reader is re- 
ferred to Table I of Ref. 22), such that the Luttinger 
coefficients all closely approach 1. Eq. 22 is then the 
sum of three identical terms. In particular, one finds 
that Im[M(u;)] cx U 2 uj 2K ^- 1 T 2 (-K p+ ) sin 2 (7rif p+ ). In 
the opposite limit, the strong coupling SO(8) fixed point 
corresponds to < g w c 

^ u 1133 ~ u 1331 ~ " C -- 

remains close to 1, spin and charge gaps drive K p+ ,-^— , 
K CT ± towards 0. If this limit were strictly attainable, 
then (73 = 0, and Eq. 22 again would become the sum 
of three identical terms, with a naive uj~ 3 dependence to 
Im[M(tj)], although it is essential to note both that: the 
coefficient of such a term would also be highly frequency 
dependent, and for the small frequencies required for the 
validity of our theory, M(w) would outgrow the pertur- 
bative regime required for the memory function's linear 
response to be valid (despite the smallness of G) a little 
before this point. 

Unlike at the strong coupling point, in the high fre- 
quency limit, 3-band umklapp processes are also relevant, 
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FIG. 5: Using the memory function, we can extract the fre- 
quency dependence of the "high frequency" optical conductiv- 
ity (in S.I. units). Here we have removed the — dependence for 
U/t ~ 0.011. Notice that one still has substantial frequency 
dependence in the crossover region and that the magnitude is 
quickly killed at high frequencies as required to satisfy the op- 
tical conductivity sum rules. The asymptotic high frequency 
dependence is o(lu) ~ oj -1 006 while as one approaches the 
Mott gap this changes dramatically, in the small frequency 
limit here o-(tj) ~ lu~ 2 ' 3 . In the inset, we show U/t ~ 0.11, 
to emphasize the small frequencies are an artifact of the weak 
coupling limit. The choice t=0.25 meV=4 * 10 14 Hz has been 
taken, such that we are restricted to oj < 2 * 10 14 Hz. 



and contribute additional terms to both Eq. 21 and Eq. 

2 



22. Letting; 



,2K, 



K„, K x -> K v and# p+ 



The real part of the optical conductivity arises as, 

5ie <r 2band (uj) = 

^e 2 u p K p ^m(M 2 band(^)) ,„„> 

H{u J + ^e{M 2band {ioW + {^m{M 2band {u J Wy [ ' 

the condition 0.2 > ui^ 1 \M 2 {, an( i(uj)\ approximates the re- 
gion where this formula can be applied. On a 3-leg ladder 
the region about which our linearization is valid is ~ | 
of the band- width (4t). This gives a high energy cut-off 
on the frequencies we can describe for the 3-leg ladder: 
oj < | so that one can describe a large range of frequen- 
cies (see Fig. 5). This constraint becomes more severe 
for N large-for N=16 one is restricted to w < 0.08t-but 
is expected to still admit a region where the memory 
function returns a(ui) cx — . 
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FIG. 6: The effects of temperature on a half- filled single chain system in weak coupling. Here, r = 27re ' /T, the 1/uj term is 
cut off at high frequencies by K p — > 1, where K p and u p are the Luttinger parameters of the single chain, and at intermediate 
temperatures we would expect to see a partial shift of weight from above the gap to form anw = Drude peak. A, B and c 
are constants, g p is the strength of the the umklapp scattering as denned in Eq.9 and 21, while ao is the weight 2 e 2 u p K p /h. 
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FIG. 7: The evolution of the two-leg ladder optical conductivity as a function of temperature, were the frequencies available 
experimentally, (top) The majority of the spectral weight would lie in a sharp peak at V3A; the onset of 2-particle scattering 
would be as \/u> — 2A; the decrease up to 3A would go as u>~ 3 , and cross over to the w _1 form of Giamarch: 19 , cut off at high 
frequencies by the larger of the resistive umklapp scattering terms, when K x — > 2, (where K x is defined in the text), such that 
the optical conductivity sum rule is satisfied. Note that this need not be a simple cross-over of exponents, one might find that 
this asymptotic behavior arrives higher than the exact low frequency result, leading to a second bump at higher frequencies 
as outlined with dashed lines. Less weight would be attributed to the exciton as the temperature approached that of the gap, 
with increasing weight at ui = 0. 
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FIG. 8: The evolution of the three-leg ladder as a function of temperature. 



